Dynamical charge susceptibility in layered cuprates: 
the influence of screened inter-site Coulomb repulsion 
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The analytical expression for dynamical charge susceptibility in layered cuprates has been derived 
in the frame of singlet-correlated band model beyond random-phase-approximation (RPA) scheme. 
Our calculations performed near optimal doping regime show that there is a peak in real part of the 
charge susceptibility \(c\,uj) at Q = (it, ir) at strong enough inter-site Coulomb repulsion. Together 
with the strong maximum in the Im x(Q, u)) at 15 meV it confirms the formation of low-energetic 
plasmons or charge fluctuations. This provides a jsutification that these excitations are important 
and together with a spin flcutuations can contribute to the Cooper pairing in layered cuprates. 
Analysing the charge susceptibilitiy with respect to an instability we obtain a new plasmon branch, 
Wq, along the Brillouin Zone. In particular, we have found that it goes to zero near Qcdw ~ 
(27r/3,27r/3). 



I. INTRODUCTION 

The phenomenon of high-temperature superconductiv- 
ity can not be completely understood without clarifica- 
tion of 'pseudogap' features seen by many experimental 
techniques in the normal state of underdoped cuprates. 
There are few informative reviews and papers devoted to 
this problem( see for example One of the most rea- 

sonable point of view is that a pseudogap in the density of 
states shows up due to an instability in two-dimensional 
copper-oxygen planes. The important question to answer 
is what is the physical origin of this possible instability 
in layered cuprates. In a series of works Di Castro and 
co-workers (see for example ||) in a frame of the phe- 
nomenological quantum critical point (QCP) hypotethis 
have suggested to explain the pseudogap as a result of 
the formation of the incommensurate charge density wave 
(ICDW). In recent paper the pseudogap phenomena has 
been considered in a context of a new hidden order pa- 
rameter with orbital antiferromagnetism (special case of 
d-wave wave order parameter) || . Constructing the in- 
tegral equation for the extended charge density waves 
instability in a frame t-J Hamiltonian proposed by An- 
derson it is easy to prove that for q = (tt,tt) the or- 
der parameter for charge instability indeed has a d-wave 
symmetry. However, the order parameter is imaginary 
H and, therefore, is not a usual charge density waves. 
This reminds the old problem of the itinerant currents 
studied many years ago in semiconductors Jj|[lO| . 

Perhaps, the most convincing evidence about the insta- 
bility in layered cuprates can be obtained directly within 
the analysis of the dynamical charge and spin suscepti- 
bilities. In contrast to the spin susceptibility the role 
played by the instabilities in the charge channel with re- 



spect to the pairing and pseudogap formation are not 
clear. Moreover, the form of the charge susceptibility be- 
yond RPA scheme was much less investigated. In partic- 
ular, there are no agreement about its analytical form if 
the effect of strong electron correlations is taken into ac- 
count (for example, compare the results in Refs. [pd|-p^[. 
Therefore, the theoretical investigation in this direction 
becomes very actual. 

Here, using the singlet-correlated band model ( 



14 



15| 



and references therein) which under simplified assump- 
tion about energy dispersion and screened Coulomb re- 
pulsion is equivalent to the t — J model |H| using the 
projecting Hubbard-like operators [[l7| we derive the ex- 
pression for the dynamical charge susceptibility beyond 
RPA-scheme. We analyze the obtained expression nu- 
merically and find that the charge susceptibility shows 
a peak at Q =(tt, 7t) similarly to a spin susceptibility 
only if one takes into account the strong enough inter- 
site Coulomb repulsion. Investigating the denominator of 
the charge susceptibility we have found the possible dis- 
persion curve of CDW-like excitations u q , for optimally 
doped cuprates. 



II. MODEL HAMILTONIAN AND 
BACKGROUND 

The Hamiltonian in the singlet-correlated band model 
in terms of Hubbard-type projecting operators, = 
\i, a >< i, 0\ reads as 



h = J2t ij *f' a *? dp + £ Mists,) - 2^] 



(1) 
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where ty is a hopping integral and Jy is the superex- 
change constant of the copper spins, a = ±1/2. Symbol 
pci-corresponds to a Zhang-Rice singlet where one hole 
is placed on copper and second distributed on neighbor- 
ing oxygen sites JlJJ] and carrying spin states, respec- 
tively, gij is a screened inter-site Coulomb repulsion of 
the doped holes. 1 + = £ H>f* + 2^ d ' dp is equation 
that determines the chemical potential. The W operators 
obey the specific commutation relation and on-site mul- 
tiplication rules of Hubbard-like operators which can be 
found somewhere [|l8],[l9) . The spin and density operators 
are expressed by projecting operators as follows 



A(k,q,w) = 



B(k,q,w) 



1 



w§ - e k+q + e k 
1 



^q - £k+q + £fe 

and frequencies are determined by the relations 



(6) 



tk+qllk+q — ,Jq n ™fc — ™fc+g 

2 ^ w c „ - e k+ a + e k M 9q) ^ 



1 = 



E 



tk+qllk+q — tkUk 



U 'q _ £ fc+9 + € k 

(7) 



S+ = S 



_ 1 T "fc- 



(8) 



J. _ ^P<Mp _ ^0,0 



where = ('J 



is a partion number, = P p dtk 



Let us now make some preliminary discussion about the 
nature of density wave formation in a frame of our model. 
We are going to study the objects which can be described 
by the following operators 



(2) is an energy dispersion of the itinerant carriers, P p d — 
is Hubbard-type bandwidth narrowing factor. The 
Fourier transform of the integrals tfc, J qi and g q are de- 
termined by the usual expressions 



i5>(k,q,c,)[^4* 



k 



1^ = 2ti (cos k x + cos ky) + 4i 2 cos k x cos k y + ■ 
J q = 2 Ji (cos q x + cos q y ) + 4J 2 cos q x cos q y + 
^l+q^k P % ( 3 ) G « = 2Gi(cosfe + cosq y ) + 4G 2 cos ft, cosq y + 



(9) 
(10) 

(11) 



6, = i£i?(k,q,a/p 



(4) 



where A(k, q, ui) and _B(k, q, are charge- and spin- 
excitation amplitudes, respectively. The frequencies of 
the harmonic motion of these objects can be found in a 
usual manner 



(5) 



Real part of 7y g corresponds to usual charge density waves 
(CDW) order parameter whereas its imaginary part is re- 
sponsible for the charge current (JC) formation pO| , pO| . 
The latter also looks very similar to the proposed earlier 
flux phase plj and orbital antiferromagnetism ||. Real 
part of £,j represents the well-known spin density waves 
(SDW) and its imaginary quantity closely related to the 
so-called spin current (JS) formation fici| , ^0|| . Symmetry 
aspects of these interesting problem has been discussed 
very recently by Nayak |22| and therefore we do not touch 
this problem here. One could immediately see from (3) 
and (4) that the equations for the frequency determi- 
nation of the real and imaginary parts looks very simi- 
lar and seem to be a good starting point to study both 
the collective charge and the spin excitations in layered 
cuprates. 

Calculating the commutators using the simplest Hub- 
bardl decoupling scheme [jlTj] (for details see the next 
section) one can find that A(k, q,w) and B(k, q,w) can 
be taken in a usual form: 



where ti(Ji, g\), i 2 (J 2 , #2), and £3 refer to the first, sec- 
ond, and third neighbors respectively. For simplicity we 
put the lattice constant equals unity. The choice of pa- 
rameters for the hopping integrals is determined by the 
right Fermi surface topology and presence of Van-Hove 
singularity in the vicinity of the Fermi level as observed 
in the experiment p3]. This satisfies at t\ = 72 meV, t 2 
= 0, and t3 = 12 meV. The screened Coulomb repulsion 
of the doped holes and superexchange integral between 
nearest copper spins were taken as G\ = 70 meV and 
J 1 = 135 meV in agreement with first principles calcu- 
lation [^if and inelastic neutron scattering experiments 
p5| respectively. 

Regarding Eqs. (0) and (||) we would like to remark 
two important features. Both equations determine the 
conditions of the instability of the system with respect 
to mentioned types of charge and spin instabilities. In 
other words Eqs. (Q) and (||) determine the divergences 
of charge and spin susceptibilities in the normal state. 
In contrast to the usual RPA-type of expression the first 
terms in (0) and (|J) proportional to the hopping inte- 
gral result from the strong electron correlation effects 
due to no double on-site occupancy constraint [|19[. It 
reflects the existence of very large on-site Coulomb repul- 
sion leading to a renormalization and enhancement of the 
susceptibilities (or in other words the presence of strong 
electron correlation). This non- Fermi liquid correction 
was firstly deduced by Hubbard and Jain |^6| for the 
spin susceptibility. The second remark is that as one can 
see from Eqs. (^) and (|^) the screened inter-site Coulomb 
repulsion plays no role for the determining of spin disper- 
sion in contrast to charge one. Therefore it becomes clear 
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that the equations for charge and spin collective excita- 
tion are quite different in presence of inter-site Coulomb 
interaction than it was widely implied. However, our 
present approach is quite preliminary and in order to see 
the difference between charge and spin collective excita- 
tions more in details we will calculate dynamical charge 
susceptibility and compare its form with spin counterpart 
obtained earlier |27|1. 



III. DYNAMICAL CHARGE SUSCEPTIBILITY. 
IMPROVED DECOUPLING SCHEME FOR THE 
EQUATION OF MOTIONS 

We start from the definition of the Fourier transform 
of the amplitude of the density operator 



&kq 
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where the Fourier transform defined as 



(13) 



In the next step we derive the equation of motion for the 
two retarded Green's functions 



(14) 

and 

i 

u c q ((ek q \e- q ))u> = —[e kq ,e- q ] + {{[e kq H]\e^ q )) u . (15) 

We calculate the commutators using the on-site repre- 
sentation of the operators. For example for the kinetic 
part of the Hamiltonian H t , it has a form 

i 

- Y.tiAyf' 1 ^ + *f d ' T (*l' T + ^f'P d )}^P d . (16) 
i 

This exact result leads, however, to the appearance of 
Green's function of higher order. Therefore we made the 
decoupling in the few steps. First, using the equation for 
the determining of the chemical potential and definitions 
of the density and spin operators determined in previous 
section one can find 



$ n , ypdpd = 1+i , fj , s z 



(17) 
(18) 



where S is an average number of the doped holes per one 
unit cell and ei its modulation. The next step is to differ 
the cases i = j and i ^ j. In latter the right hand side of 
( |l6| ) can be approximated within the usual approxima- 
tion technique 



^,{*}'WM 



1,pd\ 



1 + 5 



pd,] ,T,pd,] 



i 

+ Ltl^Pddypdd + + s *){$f d ^* d )}. (19) 
2 J 2 J 

where sf is zero in the paramagnetic phase. 

As for the case i = j it can be calculated exactly as 

[yf'P d ,H t } ^J2 t ^f ,1 ^l pd + y pd ' L ^i pd ) 

- X>(*f ,T *l ,pd + ,l ^i ,pd ) (2°) 



Using the Fourier transformation ( p"3| ) we can write the 
commutator ( jlrj ) in the equation of motion as follows 

[ekq, H t ] = Pd(tk - tk+ q )ek q 
7:7 J2[tk+ g (nl +q + n[ +q ) - t k {n\ + n l k )]e q 



AN 
2 - P, 



V 



(21) 



The last term in (|T]) still leads to the high order Green's 
function but hopefully it can be reduced to the old one 
with the help of (M), which yields 



f((e q \e- q )) = - t k+q )((e k , q \e_ q )). (22) 



Turning to the superexchange and Coulomb screening 
parts of the Hamiltonian, Hj and H c , we notice that op- 
erator \frP d,pd commutes with these terms of the Hamil- 
tonian and therefore i — j case docs not contribute to 
the equation of motion. This is in contrast to the spin 
susceptibility where the additional non-commutation of 
spin operator 5+ led to the fact that ordinary Hubbardl 
decoupling scheme is not enough and additional consid- 
eration is required Strictly speaking it means that 
although spin and charge degrees of freedom are still cou- 
pled as in the case of the usual Fermi-liquid, however, 
the simple relations of RPA-approximation of the usual 
Fermi liquid approach do not hold. 

For the case of i ^= j in site representation we have 
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- J j7 {^f* T ^j pd *P - ^ d %j pd *W}. (23) 



whereas for the case i = j 



(24) 



Restricting ourself an absence of long range spin order 
we can wright 

[e kg , Hj] = — J q {[n\ + n[] - [n\ +q + n l k+q ]}e q . (25) 

For the Coulomb screening part of the Hamiltonian the 
result reads 

[ekg, H c ] = -^g q {[n\ + n{] - [n\ +q + n{ +q ]}e q . (26) 



Substituting (El|,E5j,E6j) into the equation of motion (|2 
and taking into account the absence of the external mag- 
netic field we get finally 

Uq((ekq\e- q )) = — {n k+q - n k ) + P pd (tk - t k+q ((e k q\e-q)) 
Ztt 

+ [^(tk+ q n k +q - t k n k ) + (Gq - -^){n k+q - n k ) 



(2 - P pd ) 



-Uq]((eq\e-q)) 



2 

where n k = \{n^ + nj,). Performing a sum over k we 



(27) 



IV. NUMERICAL RESULTS AND DISCUSSION 

Let us first analyze the behavior of the obtained ex- 
pression for the dynamical charge susceptibility (pq). In 
Fig. 1 we present the calculated Re Xch(<l, 0) for the set 
of parameters described above. The peak around {tt, tt) 
reflects the nesting properties of the Fermi surface en- 
hanced by the RPA-type of the denominator. The struc- 
ture of real part of the charge susceptibility at zero fre- 
quency looks very similar to the spin counterpart. Orig- 
inally it comes from the nesting properties of the Fermi 
surface resembles by xo- However, the additional en- 
hancement due to the denominator is quite different ori- 
gin in both cases. In the spin susceptibility it results due 
to superexchange interaction having a maximum at (tt, 
tt). In the charge susceptibility the contribution from 
J q has a different sign and therefore should suppress 
the peak. On the other hand inclusion of the inter-site 
Coulomb screening repairs the situation and again leads 
to the strong commensurate peak around (tt, tt). As we 
will see later this peak is strongly dependent only on 
the value of inter-site Coulomb repulsion. This is quite 
remarkable fact if one remembers that for example the 
calculations of the charge and spin susceptibility in the 
frame of one-band Hubbard Hamiltonian with on-site U 
or ordinary t — J model show that charge susceptibility 
is 



obtain the expression for the dynamical charge suscepti- 
bility in a form 

X«Mq,wj l+{Gq _ i) Xo (q,u/) + X i(q,w) - (1 - 

(28) 

where \i (li w ) is an ordinary Pauli-Lindhard response 
function 



(q,w) 



Xo (q, w ) = -]T 



n k — n k+q 



N^uj + i0+-e k + e k+q 



(29) 



The most interesting \i (q> w ) ancl ^(qi w ) terms describe 
the contributions due to no-double occupancy constraint 
and projecting nature of Hubbard operators: 

1 tfcnfc — t k+q n k+q 

2N ■ 




(0,0) 



(0,i() 



(JMU) 



(0.0) 



xi(q,w) 



ui + i0+ - e k + e k+q 



N ^ lj + i0+ - e k + e k+q 



(30) 



1 



E 



(31) 



The comparison of the obtained equation with our pre- 
liminary consideration resulting in Eq. (Q) shows that 
important ^ 2 ^ pd ^ z(c\, lu) term has been appeared on the 



right hand side of (|31|). To our knowledge this term was 
not pointed out for the charge susceptibility so far. It 
acts as a frequency dependent "molecular field" caused 
by the transfer hopping term of the Hamiltonian (Q) . 



FIG. 1. Momentum dependence of the Re x(q, 0) through the 
route of the Brillouin zone (0, 0) — (n, 0) — (tt, it) — (0, 0) at 
T=100 K. 



rather small at Q—(tt, tt) in contrast to the spin suscep- 
tibility ]2^] . This often led to the exclusion of the charge 
degrees of freedom from the pairing interaction. How- 
ever, as it is seen here the significant contribution to the 
charge susceptibility comes from the inter-site screened 
Coulomb repulsion between doped holes. Since this in- 
teraction plays the most important role in underdoped 
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case one would expect the significant contribution to the 
pairing interaction from the charge susceptibility too. 

This fact is also seen from Fig. 2 where we present the 
imaginary part of the charge susceptibility at Q=(7r,7r) 
with and without screened Coulomb repulsion G\ . With- 
out inter-site Coulomb repulsion one sees that charge sus- 
ceptibility shows no features at low energies. This indi- 
cates that superexchange mechanism itself cannot lead 
to any charge instabilities in the system. The situa- 
tions changes drastically if one switches on the screened 
Coulomb repulsion between doped holes. At G\ = 70 
meV charge susceptibility shows a strong peak at oj c f = 
15 meV. This is quite comparable with typical spin fluc- 
tuation frequency, u s f that results from the approximate 
position of the 



J, 



1 + (G q - ^)xo(q, «) + Xi(q, w) - (1 - ^)*(q, «) = 

(32) 

through the route of the Brillouin zone (0, 0) — (it, 0) — 
(-7T, 7r) — (0,0). One has to notice, however, that the so- 
lution of Eq. (|3^) can only be obtained at larger than 
70 meV values of inter-site Coulomb repulsion. In par- 
ticular, at the optimal doping, i.e. p = 0.16, the stable 
solution in the whole Brilloune zone was found only at 
Gi = 150 meV. In Fig. 3 we present the resulting en- 
ergy dispersion of CDW frequency, uj^ dw at tempera- 
ture T=100 K. As one can see there are three discon- 
nected branches of the frequency curves. It happens due 
to the fact that instability occurs not along the whole 
Brllouin zone but rather at some critical wave vectors. 
In particular, it is clearly seen that the charge system 
is very unstable near the incommensurate wave vector 
(2/7r/3, 27r/3) where curve is crossing zero. This means 
that the possible CDW must occur at this incommensu- 
rate wave vector Qcdw ~ (27r/3, 2n/3). This quite a 
remarkable result confirms the difference between para- 
magnon charge fluctuations which are commensurate and 
real CDW order parameter having incommensurate wave 
vector. 
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FIG. 2. Frequency dependence of the Im x(Q> w) at Q = 
(7T, tt) for Gi=0 (long-dashed curve) and Gi= 70 meV (solid 
curve) at T = 100 K. The arrow indicates the approximate 
position of the plasmon or charge fluctuation frequency w c / 
ps 15meV. 



peak in Im Xs, ^sf and which is also in order of 15 
meV in layered cuprates [M. On the basis of Figs. 1 
and 2 we can conclude that charge susceptibility may 
play an important role in the scattering mechanism and 
in the Cooper pairing instability correspondingly. How- 
ever, such a large value of inter-site Coulomb interaction 
one could expect only at low doping level. Therefore, 
we think that the spin fluctuation scenario of supercon- 
ductivity is much more robust than charge fluctuation 
one in high-T c cuprates. One would expect, however, a 
strong charge density wave like instability in underdoped 
cuprates since inter-site Coulomb repulsion is the most 
important at low doping concentration. 

In order to analyze these possible instabilities in the 
charge subsystem we also investigate the charge suscep- 
tibility with respect to a CDW formation. The CDW 
frequencies range can be obtained analyzing the denomi- 
nator of dynamical charge susceptibility expression ( 28 ) . 
Therefore, we solve the equation 




(0,0) 



(0,0) 



FIG. 3. Momentum dependence of the CDW frequency acous- 
tic mode at Gi= 150 meV. 



V. CONCLUSION 

In a conclusion we have obtained in a frame of singlet- 
correlated band model beyond RPA approximation the 
new analytical expression for the dynamical charge suc- 
septibility. We have shown that it produce a strong peak 
at Q = (tt, 7t) at large enough screened Coulomb re- 
pulsion between doped holes. This may result in a sig- 
nificant contribution to the pairing interaction from the 
charge susceptibility at low doping level where Coulomb 
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repulsion plays an important role. The analysis of the 
instability of the system with respect to a CDW forma- 
tion shows that the system is the most unstable at the 
incommensurate wave vector Qcdw ~ (27r/3, 2tt/3). 
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